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Abstract: We provide the formalism for the quantization of systems of coupled 
bosonic and fermionic fields in a time dependent classical background. The occupa- 
tion numbers of the particle eigenstates can be clearly defined and computed, through 
a generalization of the standard procedure valid for a single field in which Bogolyubov 
coefficients are employed. We apply our formalism to the problem of non-thermal 
gravitino production in a two-fields model where supersymmetry is broken gravi- 
tationally in the vacuum. Our explicit calculations show that this production is 
strongly suppressed in the model considered, due to the weak coupling between the 
sector which drives inflation and the one responsible for supersymmetry breakdown. 
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1. Introduction 

The analysis of quantized systems in a classical background can be very useful for the 
study of various phenomena that arise in quantum theories, as for example particle 
production. The study of matter in external electromagnetic fields dates back to 
the first years of quantum field theory P, |3| . For what concerns gravity M , the semi- 
classical approximation is often compulsory, due to the lack of a consistent quantum 
theory. Despite of this, it turned out very successful in describing phenomena as 
particle creation from black holes |5| or the generation of the perturbations in the 
inflationary Universe ||. 

In the last ten years, this semiclassical approach has been applied to non-thermal 
particle production after inflation. In this case, the classical background is given by 
the inflaton field, which is coherently oscillating about the minimum of its potential. 
The first analyses of this phenomenon were performed in J7|, but its full relevance 
was appreciated only a few years later in the case of production of scalars p, |9], |10| ] . 
In the work M this non-perturbative production has been called "preheating" , since 



it is usually followed by a ordinary phase of perturbative reheating. It has been 
there understood that preheating of bosons is characterized by a very efficient and 
explosive creation, even when single particle decay is kinematically forbidden. This is 
due to the coherent inflaton oscillations, which allow stimulated particle production 
into energy bands with very large occupation numbers. 

Less attention was initially paid to non-perturbative production of fermions. 
Indeed, the efficiency of this process seems to be strongly limited by Pauli blocking, 
which does not allow for occupation numbers bigger than one. However, also this 
production turned out very relevant, as the first complete calculation [11] of the 
inflaton decay into heavy (spin 1/2) fermions during preheating showed. 1 Indeed, 
if one only considers the most natural interactions ip ip and </> 2 \ 2 of the inflaton 
cf) to fermions ip or to bosons Xi fermionic production occurs in a mass range much 
broader than the one for heavy bosons, and this can "compensate" the limit imposed 



by Pauli blocking |TT 



These results were soon applied fT4] , 15] to non-thermal gravitino production, 
since the equations for the different components of the gravitino field can be reduced 
to the one of a spin 1/2 particle. As had also been realized in [IB|, the transverse 
gravitino component is always very weakly coupled to the background and decoupled 
from the other fermions, so that the production of its quanta is negligible. However, 
the works ]14|, [L5j also studied the production of the longitudinal component, con- 
cluding that it easily exceeds the limits imposed by primordial nucleosynthesis (the 
so called "gravitino problem" fll7|). The analyses of |L4], [L5| were extended in ][L8[ [19]] 
and followed by several related works [^, |2l| . 

Most of these analyses support the conclusions of []14 , |i"5f of a gravitino over- 
production at preheating. However, in all the explicit calculations, only the case of 
one chiral superfield with supersymmetry unbroken in the vacuum was considered. 
This last issue is however crucial to understand the gravitino production: in the 
super-higgs mechanism (in the unitary gauge), the gravitino longitudinal component 
is provided by the goldstino, which is present only when supersymmetry is broken. 
As pointed out in the works [B|, during the cosmological evolution supersymme- 
try is broken both by the kinetic and the potential energies of the scalar fields of the 
theory. However, these fields are expected to be settled in their minima now. If in 
these minima supersymmetry is unbroken the gravitino has only the transverse com- 
ponent and it is massless. The calculations performed in these schemes show that 
one fermionic component is produced at preheating. However, we have remarked 
that at late times it cannot be the longitudinal gravitino component. Since this 
fermion is the partner of the inflaton field (in case of only one chiral multiplet the 
scalar is necessarily the inflaton), it should be better denoted as "inflatino". Since 



1 Among other interesting studies on production of fermions (not all of them related to preheat- 
ing) we mention [ff^, E3. 



this field does not have a gravitational decay rate, 2 we conclude that preheating does 
not contribute to the gravitino problem in models with supersymmetry unbroken in 
the vacuum. 

One is thus led to consider schemes where both the issues of inflation at early 
times and supersymmetry breaking today are included. From COBE normalization 
of the scalar metric perturbation, the relevant scale of inflation (when only one field 
is present) is expected to be about 10 13 GeV. Supersymmetry provides instead a 
solution to the hierarchy problem if it is broken at about the TeV scale. Although 
in principle one may construct a model where a unique field satisfies both these re- 
quirements, we do not consider this option as the most natural one. Moving to the 
two fields case, the simplest possibility is to consider two well separate sectors, one 
of which drives inflation, while the second is responsible for supersymmetry breaking 
today. We have in mind a situation in which no direct coupling is present between 
the two fields in the superpotential, so to have a strong suppression in the inter- 
actions between the inflaton and the field which provides the longitudinal gravitino 
component. As we describe below, with two chiral supermultiplets the longitudinal 
gravitino component is coupled to one other fermionic field (the matter component 
orthogonal to the goldstino). For simplicity, in the rest of the introduction we refer 
to these two fields simply as to the gravitino and to the inflatino, although this in- 
terpretation is true only in the vacuum of the theory and more rigorous definitions 
will be provided below. The coupled system is particularly involved, so that trying 
to guess its behavior without an explicit calculation results very difficult. One can 
guess that at the preheating era only the inflatino field is produced. However, there 
is the potential worry that much later (on a physical time-scale of the order the 
inverse gravitino mass, when supersymmetry is equally broken by the two sectors of 
the theory) a fraction of inflatinos is "converted" in gravitinos. As remarked in [|T^] . 
this worry requires an explicit calculation in a full supergravity context. 

The problem of gravitino production is thus reduced to the problem of a (quite 
involved) coupled system in the external background constituted by the scalar fields 
of the theory. The most difficult part of this analysis is to provide a formalism in 
which the coupled system is quantized, with a clear definition of the occupation 
numbers for the physical eigenstates. This is a very interesting problem in itself, 
which can have several other applications besides the one we will consider in the 



present paper. In the one field case, the procedure is well established |12], p2| . One 
first quantizes the system and expands the canonical hamiltonian in the creation 
and annihilation operators of the field. The evolution of the background creates 



2 As also remarked in | pT| , supersymmetry requires the inflatino decay rate to be comparable 
with the inflaton one, such that the inflatino is expected to decay at reheating (this is actually 
true provided that supersymmetry breaking at that time is sufficiently small). In ref. EJ , explicit 
calculations are however performed only with one relevant superfield and unbroken supersymmetry 
in the vacuum. 



a mixing between the positive and negative energy solutions of the field equation, 
which has the consequence of driving the hamiltonian non diagonal, even if one 
takes it to be diagonal at initial time. A diagonal form is achieved through a (time 
dependent) redefinition of the creation/annihilation of the fields. The two coefficients 
of this diagonalization are known as Bogolyubov coefficients and can be easily related 
to the occupation number for the quantized field (consistency of the quantization 
requires a relation between these coefficient; this condition is however preserved by 
the equations of motion of the system). 

In the first part of the present work we generalize this procedure for systems 
of multi- fields, both in the bosonic and in the fermionic case. Although far from 
trivial, this generalization can be presented in a remarkably simple form. By choosing 
a suitable expansion of the fields we can repeat each step of the above analysis 
substituting the Bogolyubov coefficients with two matrices a and (3 (the expansion 
is now performed in a basis of creation/annihilation operators, each corresponding 
to a physical eigenstate of the system after the diagonalization of the hamiltonian). 
We can obtain a system of first order differential equations for these matrices; the 
condition on a and (3 for a consistent quantization are also very simple and are 
preserved by these equations. Finally, the expression for the occupation numbers 
is an easy generalization of the one valid in the one field case. 3 This first part is 
divided in two sections, the first of which is devoted to bosons, while the second 
one to fermions. This second section is further divided in two parts. In the first 
one we consider the case in which the fermionic fields are coupled only through the 
"mass matrix" , while in the second one we consider a more general system which is 
necessary for the application to the gravitino case. 

The second part of the work is completely devoted to this application. In sec- 
tion |4.1| we introduce the quantities relevant for the calculation, following the nota- 



tion of [ 19| . In section fO| we describe the model we are considering, discussing the 



evolution of the scalar fields. In section [4.3| we write the effective lagrangian for the 
longitudinal gravitino component and the matter fermion in the two chiral supermul- 
tiplets case. In the last two subsections we present our results for the gravitino and 
the inflatino production. In section |4|4] we present analytical results for the case in 
which supersymmetry is unbroken in the vacuum, while numerical results for the case 
of broken supersymmetry are shown in section |4.5| . As we will see, the final gravitino 
production is strongly sensitive to the size of the supersymmetry breaking: our nu- 
merical results indicate that this production goes to zero in the limit of a vanishing 
final supersymmetry breaking. This limit is in agreement with the analytical result 
of section [4.4| and with the fact that the longitudinal gravitino component is actually 
absent when supersymmetry is preserved. We thus conclude that in the model con- 



3 As in the one field case, a state of vanishing initial occupation number can be easily defined, 
provided the system evolves adiabatically at initial time. 



sidered here, non-thermal gravitino production is strongly suppressed. This appears 
as the consequence of the weak coupling between the two sectors responsible for in- 
flation and supersymmetry breakdown, and of the strong hierarchy between the two 
scales that characterize them. 

A short description of our results for the gravitino and inflatino production can 
be found in J2|. 



2. System of coupled bosonic fields 

In this section we consider the coupled system of N bosonic fields {0j} in a FRW 
background described by the action 



1 



S = - I d x d^i d^(j)i - m\ , fa + £ R . 



XV, 



We use conformal time rj, such that the metric and the Ricci scalar are g^ v = 
a 2 {rj) diag(l, —1, —1, —1) and R = —6 d/a 3 , where a is the scale factor of the Uni- 
verse and dot denotes derivative with respect to rj (summation over repeated indices 
is understood). The last term describes a possible non- minimal coupling (£ ^ 0) of 
the scalar fields to gravity. 

The (symmetric) mass matrix mfj is assumed to be a function of some external 
(background) fields. The only assumption that we do on these external fields is 
that they are constant (or better, adiabatically evolving) at the very beginning 4 and 
at the very end of the evolution of the system. In these regimes, the matrix m 2 j 
becomes also constant and the fields which diagonalize it become free fields, whose 
masses are precisely given by the eigenvalues of mfj. However, during the evolution 
the different entries of m 2 - are allowed to vary, and the (time dependent) eigenvalues 
of mfj are interacting fields whose masses change in time. These masses can change 
non adiabatically and this may in general lead to particle production. The aim of 
this section is to give a precise definition of the occupation number and to provide 
the formalism to calculate it. 

It is most convenient to consider the "comoving" fields </?j = a (pi. For these 
fields, the above action (12.11) rewrites 5 



S = ^ J d 4 x [<pi cpi - pi fi?- (fj] 



fi* = a 2 ml + ( -A + - a (6£ - 1) ) , (2.2) 



We require an initial stage of adiabatic evolution to consistently define vanishing occupation 
numbers for the bosons at initial time. 

5 We do not necessarily need a cosmological motivation for the analysis that we perform in the 



rest of this section. The action (2.2) could indeed also arise in flat space, with a non-diagonal flfj 
coming from some general interactions between the bosons <fi and some other fields which have 
been integrated out. 



where A is the laplacian operator. We can also write the hamiltonian of the system, 
which, in terms of the fields ipi and their conjugate momenta 

dC 



n, 



d(fi 



reads 



(2.3) 



[2A) 



The frequency matrix Vl 2 - which enters in the above expressions is also in general 
time dependent and non diagonal. At any given time, it can be diagonalized with 
an orthogonal matrix C 



C T {rj) tt 2 {rj) C {rj) = to 2 {rj) diagonal . 



(2.5) 



We denote by (p = C T <p the bosonic fields in the basis in which the frequency matrix 
is diagonal. We also denote by uf the i-th entry of the diagonal matrix uj 2 . The set 
of uJi represents the energy densities of the (time dependent) physical eigenstates of 
the system (p^ 

We now show that the occupation numbers of these fields can be defined and 
computed generalizing the usual techniques based on Bogolyubov coefficients valid 
in the one field case. The first step to do this is to consider a basis for annihila- 
tion/creation operators {<2j} and {a\} and to perform the decompositions 



'.1 



d 3 k 



a 



(2tt 
d 3 k 



3/2 



(2tt 



3/2 



i kx 



i kx 



hjk (v) a k (k) + e ik -*h*(r))ai(k) 



h jk tn) a k {k)+e- ik - x h* jk {rj) a\{k) 



;2.6) 



The reason why we explicitly factorized the matrix C in these decompositions will 
be clear soon. Due to the fact that the fields are coupled together, the matrices h 
and h are generically expected to be non diagonal. 
To quantize the system, we impose 



[<Pi (x) , U 3 {y)\ = i 5 3 (x - y) 6. 



1.1 



for the conjugate fields, and 



ai (k),at(p) =5 3 (k-p)5. 



1:1 



[2.7) 



X8) 



for the annihilation/creation operators. We can satisfy both these relations requiring 



hti - h* h T 



i S 



;2.9) 



as it can be easily checked from the decomposition (|2.6). 
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From the action (|2.2|) , one deduces the second order equations of motion 



<pi + n I \ <pj = . 



(2.10) 



However, one can achieve a system of only first order equations by setting some 
additional relations between the conjugate fields ipi and IT. We want these relations 



to generalize the one which is usually taken in the one field case, see i.e. ||22|| . Also 



we want them to allow a rewriting of the hamiltonian (2.4) in a simple and readable 
form. The sets of fields where this generalization is most evident is given by {ipi, IT = 
(C T n)j}. These fields are decomposed as in eqs. (pUD, only without the C matrix 



before the integrals. In terms of these fields, the hamiltonian (^4) rewrites 



H = / G? 3 X - ( fli ft; + Ujf (pi 0i 



(2.11) 



since the frequency u is diagonal. One is thus led to impose the conditions 6 

p— i J v ujdri' i f v ujdr)' 

h = = A+ — B , 

h= ^= A + — — B (2.13) 



V 2 a; V2cj 
which are indeed a natural generalization of the one which is usually taken in the 



one field case f22|| . For one field, A and B are numbers, known as Bogolyubov 
coefficients. In our case they are N x N matrices. The analysis of the system is in 
our case simplified if we consider, rather then the matrices A and B, the combinations 



a = e -iS v »*tf A, 

= e » pudrf g _ 



(2.14) 



The above condition 



is satisfied if the matrices a and (3 obey the relations 
1. 



aa^ -(3* f 
a /3 f - (3* a T 



0. 



(2.15) 



These relations can be imposed at the initial time, and are preserved by the evolution, 
as we shortly discuss. In the one field case, they reduce to the usual condition 



a 



\A\ 



\B\ 



1. 



6 Equations ( |2.13 ) are written in matrix notation. In general, for any function f(u>i) and any 
matrix M, we use the notation 



(/ M M) tj = f (u>i) M K 



(M/(a;)V.=M ij /(^). 



(2.12) 



As we have said, the evolution of the system can be described by two sets of first 
order differential equations. The first set is obtained by inserting eqs. ( |2.6| ) into the 
definition of the conjugate momenta, eq. Q2.3Q 



h = h-Th, (2.16) 

where we have defined the matrix 

r = c T c, r T = -r. (2.17) 

The second set of equations is obtained by rewriting eqs. ( p,10|) in terms of tpi and IT 



h = -Th-uj 2 h. (2.18) 



We can now use relations ( |2.13| ) and decouple the terms proportional to a and 



(3, so to arrive to the final result 



UJ 

a = —iua-\ (3 — / a — J (3 , 

2uj 

(3 = ^a + itoP - Ja- I (3 , (2.19) 
2 ui 



where we have defined the matrices 



2 V Juj Juj 



J = U y fcT^=-^=T y fc), J T = J. (2.20) 
2 V a/w Juj 



In the one field case, I = J = Y — 0, and the above system reduces to the 
equations for the two Bogolyubov coefficients 



A = —e 2i S v " d v'B, B = —e- 2 ^^ d ^'A, (2.21) 
2uj 2uj 



already discussed in the previous literature (see i.e. ]T0[)- In the one field case the 
only source of nonadiabaticity is related to a rapid change of the only frequency uj(rj), 
so that the system is said to evolve adiabatically whenever the condition uj uj 2 
is fulfilled. In the present case, there are more sources of nonadiabaticity, related 
to the fact that now the frequency is a N x N matrix. This is associated with 
the presence of non-vanishing matrices I and J in the equations of motion for the 
matrices a and j3. 



It is a straightforward exercise to show that the above equations (|2.19|) preserve 



the normalization conditions (|2.15|) , due to the properties I T = —I and J T = J. 



In the one field case, the number of particles is given by the modulus square of 
the second Bogolyubov coefficient, \B\ 2 . We now show that also in the multi-field 
case it is generally related to the matrix (3. To see this, we decompose also the energy 



S 



density operator "H (see eq. (|2.4j) ) in the basis of annihilation and creation operators 

"=(°M(/ fl $)(:;)■ ^ 

From eqs. ( |2.6| ), one sees that the N x N matrices £ and T which enter in this 
decomposition are given by 



T = -(h J 'h + h T uo 2 h^j . (2.23) 

We can now generalize the procedure adopted in the one field case. The matrix 
that appears in eq. ( |2.22| ) can be put in diagonal form in a basis of new (time 
dependent) annihilation/creation operators. Only when the hamiltonian is diagonal, 
each pair of (redefined) operators can be associated to a physical particle, and used 
to compute the corresponding occupation number. The explicit computation gives 

E = - (a^oua + ftcup) , 

T = - (a T u (3 + (3 T u a) , (2.24) 



so that we found that expression ( |2.22| ) evaluates to 



In terms of the redefined annihilation/creation operators 7 

a) J = (p a* ) Ut 



(2.28) 



the hamiltonian is thus diagonal (remember that in eq. fl2.5|) uj was defined to be 
diagonal), and, after normal ordering, it simply reads 



H = I d 3 kuialai. (2.29) 



7 The relation (2.28) is inverted through the matrix 



at -/3t ' 
-(3 T a T ' ' 



(2.26) 



as can be easily checked from conditions ( 2.15| ). We thus see that also the relations 



a^a-P^P = l, a^j3*-fta*=0 (2.27) 

hold for the whole evolution. 



We choose at initial time a(r/o) = 1, (3(vo) = 0, so that conditions (|2.15|) are 
fulfilled. We also choose the initial state of the theory to be annihilated by the oper- 
ators dj. At any generic time, the occupation number of the i-th bosonic eigenstate 
is given by (notice that in this expression we do not sum over i) 

N l ( V ) = (a\a t ) = ((3*(3 T ) i .. (2.30) 

In the one field case the above relation reduces to the usual iV = \B\ 2 . We see 
that our choices correspond to an initial vanishing occupation number for all the 
bosonic fields. 



3. System of coupled fermionic fields 

We now consider a system of coupled fermions. We divide this analysis in two 
subsections. The first of them extends to the fermionic case the results obtained for 
bosons in the previous section. Because of the repeated analogies, the discussion is 
here shorter than the above one, where more details can be found. In the second 
subsection we study a more general system of equations, which can be also relevant 
when the background is not constant. In particular, these can be important for 
cosmology, where the expansion of the Universe provides a preferred direction in 
time. In the next section we will indeed discuss non-perturbative production of 
gravitinos as an application of this analysis. 

3.1 The simpler case 

Let us consider the coupled system of N fermionic fields {^} in a FRW background 
described by the action 



S — I (fx sf—g ipi 



V>i- (3-1) 



The gamma matrices 7 M in FRW geometry are related to those (7 M ) in flat space by 
7 M = a~ 1 7 M , where a(rj) is the scale factor of the Universe. As before, conformal 
time rj is used, and the matrix My is considered to be a function of some external 
background fields. The requirement that the action is hermitean forces M to be 
hermitean as well. For simplicity we will also take it to be real. We also require My 
to be constant (better, adiabatically evolving) at very early and late times, but we 
do not make any other assumption on its evolution. 

After the redefinitions X{ = ipi a 3 / 2 , m = aM, the action (|3.1|) reads 



3 = J d A x X t [i 5^ Y <9 M - my] X 3 , (3.2) 

leading to the equations of motion (in matrix notation) 

(iYd^-m)X = 0. (3.3) 



10 



The canonical hamiltonian is instead 



H = Jd 3 xH = J<PxX[-iY d, + m] X . 



(3.4) 



In analogy with the bosonic case, we expand the fermionic eigenstates into a 
basis of creation/annihilation operators 

-^e* k - x \u* k (k, rj) a k r (k) + Vj k (k, n) K k (-k) ' , 

(27T) 1 L J 

where the matrix C is employed into the diagonalization of the mass matrix m 



X,, (x) = C, 



(3.5) 



jj, = C T mC , C orthogonal . 
We also define the matrix (dot denotes derivative with respect to 77) 

r = c T c, r T = -r, 

and the "generalized spinors" 



(3.6) 



(3.7) 



ui j = 



vy = 



(3. 



with ip + = f M and ip- = eigenvectors of the helicity operator er • v/|v|. 

Let us consider a set of fields X { which satisfy the above equations ( |3.3|) . Due to 
the fact that the matrix rriij is real and symmetric, then also the fields Xf = C Xj 
(where C is the charge conjugation matrix 8 ) are solutions of (|3.3j ). 9 As a consequence, 
one can impose the relation U r {k) = CVj (—k), or, using eqs. (|3.8|), 



-U* 



v_ = ui. 



(3.10) 



Doing so, we have only to deal with the U± matrices. Taking the momentum k along 
the third axis, their equations of motion read 



U± = -i k C/ T =f iii U± - V U± . 
The quantization of the system requires 



a ri (k) , a\j (p) 



(3.11) 



3 In our computations, we take 



1 
-1 



a, 

~<Ti 



6^ (x 


- y) kj , 


6^ (k 


- p) 5 rs Sij 


6^ (k 


- p) 5 rs 


c = 


= ^7°7 2 =( 



ia 2 

i <T 2 



(3.12) 



(3.9) 



where a are the Pauli matrices. 

9 This may allow one to consistently define the Majorana condition Xf = Xi 



11 



We can simultaneously satisfy these conditions by setting 10 

U+Ul + Ut Ifl = 2 1 , 

U+Ul = 17117%. (3.13) 

These conditions can be imposed at initial time, and are preserved by the evolution 
of the system (as it is easily checked from eqs. ( ft.llj) ). 
We define the diagonal matrix 11 

u = ^k 2 + /i 2 , (3.14) 

and we further expand 

U + = (l + ^) 1/2 e-<r<^' A _ (l-^jP^'B 

U. = (l - £) 1/2 e-«r»*' A + (i + ^"'e'/"""'' B 

1 _^' /2 a+ f 1 + ^ I/2 /3 , W (3.15) 



so that the above conditions ( |3.13j ) are satisfied if the matrices a and (3 obey the 
relations 

a a 1 + /?* (3 T = 1 , 

aft -(3*a T = 0. (3.16) 

In the one field case, A and B are numbers, known as Bogolyubov coefficients. In 
our case they are N x N matrices. The matrices a and f3 are introduced since their 
equations of motion assume a simpler form then the corresponding ones for A and 
B. In the one field case, the above relations ( |3.16|) reduce to the usual condition 
\A\ 2 + \B\ 2 = 1. 

For fermions, the evolution equations for the matrices a and (3 can be obtained 
in a more straightforward way with respect to the bosonic case. This is because 
eqs. (|3.11 ) are already two sets of first order differential equations. Using the above 



10 Notice that all this analysis generalizes the one made in the one field case. For the latter, we 
follow pH . A detailed exposition with a notation closer to the present one is found in |M . 

11 This definition is meaningful, since both uj and [i are diagonal matrices. More simply, it can be 
understood as a relation between their eigenvalues. See also the footnote with eq. ( p. 12 ) for some 
clarification about our notation. 
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decomposition ( 3.15 ), after some algebra we arrive to the final expressions 



12 



at = [—i uj — I] a + 
fi k 



fi k 



+ J 







2uj 2 



J 



2uj 2 
a + [i uj — I] (3 



where we have defined the matrices 



21 
2J 



it \ 1/2 / it \ 1/2 

i + r i + - 

uJ V UJ J 

^1/2 

UJ, 



/i\l/2 



UJ 



r i 



it \ 1/2 / it\ 1/2 

i_^ 1/2 r(i + -^ 1/2 

UJ/ V UJ, 



(3.18) 



I T = -I 



J T = J. (3.19) 



One can easily verify that these equations preserve the above conditions ( |3.16| ). 

To properly define and compute the occupation number for the fermionic eigen- 
states, as before we expand the energy density operator TL (eq. ( |3.4|) ) in the basis of 
annihilation and creation operators 



n 



a] , bj 



Sii 
T 



-E T 



a-i 

It 



Using eqs. 



and (13.81), we find 



1 



T(r]) = - -Ul/j,U--U_iiU + + UlkU + -UtkU- 



while eqs. ( |3.15| ) lead to 



£ = oft uj a — ft uj (3 , 
T = —a T uj (3 — (3 T uj a 



(3.20) 



(3.21) 



(3.22) 



We have thus 



12 In case of only one superfield, these equations simplify to 

A = ' n. u = 



A fc _2i fudr,' R £ = A fc c -2i r^d n ' A 



2u 2 



2w 2 



(3.23) 



(3.17) 
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In terms of the redefined annihilation/creation operators 13 

(»)-(; ?)(;) (3 - 26) 

the hamiltonian is thus diagonal, and, after normal ordering, it simply reads 

H = J d 3 k u t (a\ di + l\ b^j . (3.27) 

We choose at initial time a(rjo) = l,(3(r)o) = 0, so that conditions ( |3.16|) are 
fulfilled. We also choose the vacuum state of the theory to be annihilated by the 
initial annihilation operators ai and hi. At any given time, the occupation number 
of the i-th fermionic eigenstate is given by (notice that in this expression we do not 
sum over i) 

N t (v) = (ald l ) = (blb l ) = (P*p T ) u . (3.28) 

In the one field case the above relation reduces to the usual N = \B\ 2 . We 
see that our choices correspond to an initial vanishing occupation number for all the 
fermionic fields. Notice that particles and antiparticles have the same energy and are 
produced with the same amount, due to the reality conditions that we have imposed 
on the system. Finally, we observe that the first of conditions ( |3.16| ) guarantees that 
Pauli blocking is always satisfied. 

3.2 A more general case 

We now consider a more general action for the coupled system of N fermionic fields. 
For future convenience, here we switch to the signature — , +, +, + for the Minkowski 
metric, and we then work with the gamma matrices 

in flat space. 

By a suitable conformal rescaling of the fermionic fields and of their masses as 



we did before eq. (|3.2j ), we can again remove the scale factor of the Universe from 



the kinetic term for the fermions. However, we are now interested in a more generic 



13 r 



The matrix in eq. (3.26) is unitary, so its inverse one is precisely given by 



a 



(ft 

-I3 T a T 



(3.24) 



as can be easily checked from conditions ( 3.16| ). We thus see that also the relations 



x*a + 0tp=l, a f (3* - ft a* = (3.25) 
hold for the whole evolution. 
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system, so that we consider, instead of ( ft.2| ), the action 

S = J d A x X m [ 7 °<9 + f N d, + M] mn X n , (3.30) 

where N and M are two N x N matrices of the form 

N = N 1 +fN 2 , M = M 1 + fM 2 . (3.31) 

The matrices N and M are assumed to be functions of some external fields. 
We consider a situation in which these fields evolve in time. This time dependence 
justifies the general form for the action that we want to discuss. As we will see in the 
next section, this analysis can have relevance for cosmology, where the expansion of 
the Universe provides a natural direction for time. However, the system ( |3.30 ) could 



also arise in flat space from some general interactions between the fermions Xj and 
other fields which have been integrated out. As in the previous analyses, our main 
goal is to discuss the definition of the occupation number of the physical eigenstates 
of the system, and to provide the formalism to calculate it. 

We list here our assumptions on the matrices M and N. First, we require them to 
change adiabatically at initial times, so to consistently define the initial occupation 
numbers. Then, we assume Mi — > constant, M 2 — > 0, N — > 1 at late times, so 
to recover a system of "standard" decoupled particles at the end (indeed, one can 
choose the basis of fields X% such that the matrix M is diagonal at the end). The 
requirement of an hermitean action translates into the conditions 

N\ = Ni, Mi = M\ , M 2 = -M\ . (3.32) 

For simplicity, we will only consider real matrices, so that Ni, N 2 and Mi are required 
to be symmetric, while M 2 antisymmetric. Finally, we impose an additional condi- 
tion, which is that the kinetic term for the fermions "squares" to the D'Alambertian 



operator □. If we take the equations of motion following from (|3.30| ). 

[l°d +fNdi + M]X = 0, (3.33) 
and we multiply them on the left by [7°<9 + 7* N di + M]t, we get 

j<9 2 -N ] Nd 2 i+M ] M- 7° [(doN) f di + (doM)] }x = . (3.34) 
We thus require N = 1. This rewrites on the conditions 

N? + iVf = 1 , [Nt, N 2 ] = . (3.35) 

Our strategy is to reduce this problem to the one we have already discussed. 
That is, we perform some redefinitions of the fields to put the action ( |3.30|) into the 
form ( |3.2j ), where we perform the canonical quantization of the system in the way 
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described in the previous subsection. The first of these redefinitions strongly relies 
on the above conditions (|3.35|) . If iV is a unitary matrix, we can find a hermitean 
matrix $ such that 

iV = exp (2 $ 7°) , $ f = $. (3.36) 
Due to the properties of the 7 matrix, this amounts to 

cos (2 $) = Nt , sin (2 $) = N 2 . (3.37) 

After the redefinition X = exp(— 7 $) X, the equations of motions fl3.33p acquire 
the form 

(fdo + if ki + M^X = 0, (3.38) 

where we have expanded the fermions into plane waves Xi(r], k) = e* k ' x Xi(rf) and 
introduced the new "mass matrix" 

M = exp (7 $) [M + fd ] exp (-7 $) 

= Mi + 7 M 2 • (3.39) 

Notice that the two matrices Mi and M 2 are symmetric and antisymmetric, respec- 
tively. This means that, in the one field case, one recovers the standard equation 

(fd + ifk i + m)6 = (3.40) 

for spin 1/2 fermions. 

We have instead to perform a further redefinition of the fields. 14 Setting X = L H, 
we have 

L T \l (7 do + if ki) + Mi L + 7 (do + M 2 ) L] H = . (3.42) 

The matrix L can be chosen such that (do + M<i)L = 0. In particular, since M 2 is 
antisymmetric and real, L can be taken orthogonal. The equations of motion for the 
fields 5 can be thus cast in the form 

(fdo + ifk i + L T M l l)j-~=Q, (3.43) 

that is with the identity matrix multiplying the term which depends on the momen- 
tum and without any 7 dependence in the "mass" matrix. 



14 Contrarily to naive expectations, the combination 

cos $ d cos $ + sin $ <9 sin $ C Af 2 (3-41) 

can be non vanishing at late times, even if the matrix N is approaching 1. This occurs for example 
in the application that we discuss in the next section. If in that case we canonically defined the 
hamiltonian H starting with the fields X, the term ( 3.41 ) would give H a contribution proportional 
to 7 which does not vanish at late times. The procedure described in the main text removes 
this problem. 
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These equations (and the respective action for the fields E$) are exactly of the 
form considered in the previous subsection, so that we can apply the quantization 
procedure discussed there. As before, the procedure is to canonically define the 
hamiltonian starting from the set of fields H and to expand it in a basis of cre- 
ation/annihilation operators. The occupation numbers are then calculated after the 
diagonalization of the hamiltonian. It is possible to show that this procedure can be 
carried out starting from any of the basis for the fermionic fields, once the hamilto- 
nian has been canonically defined in the basis 5. One can indeed explicitly verify 
that these calculations lead to the same results for the occupation numbers of the 
physical eigenstates, in the same way as in the bosonic case the result (|2.23| ) can be 
computed starting from each of expressions (|2.4|) or (|2.11|) . 

We thus have 



H = E 

= X 



ifh + L T Mi L 



E = X 



X 



fe .e27 * + e-7°* Ml e7 * 



X, 



(3.44) 



depending on which basis we work. In particular, when working with the X,- L or the 
Xi fields, the explicit knowledge of the matrix L is not needed. 15 We present here 
the computation in the initial basis Xi, which we found more convenient (in the 
numerical calculations) for the application that we present in the next section. In 
this basis, the hamiltonian (|3.44|) has the form 



H = X 



i f k { N + Mi + 7° M 2 



X 



(3.45) 



where the matrices M and M 2 can be obtained from eqs. (|3.44j) and ( |3.39|) . At the 
end of the evolution, we simply have Mi +7 M 2 = Mi diagonal, so that one recovers 
the "standard" hamiltonian for a system of iV decoupled fermions whose masses 
coincide with the ones of the equations of motion (which also become "standard"). 

To analyze the system during the evolution, we instead decompose the spinors 
Xi as in eqs. (|3.5|) and ( |3.8|) . 16 Taking the third coordinate along the momentum k, 
the equations of motion ( |3.33j ) rewrite: 

=F«(Mi =f iM 2 )U± -ik (N x ± iN 2 )U T . (3.47) 



U 4 



It is straightforward to check that they preserve the conditions 



u+ u\ + u*_ u 



T 



2-1 



u, w 



u*_ul, 



(3.48) 



which ensure the consistency of the canonical quantization. 



15 However it is crucial that M2 is antisymmetric, which allows L to be orthogonal. 

16 The charge conjugation matrix now reads C = — 7 7 2 , so that conditions ( 3.10| ) are replaced by 



V, = -i U*_ 



V-=iUl 



(3.46) 
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We then expand the hamiltonian formally as in eq. (|3.20| ), where now the 8 and 
T matrices read 

S = Ul k [iVi + i N 2 ] U- + Ul.k [JVi - i N 2 ] U+ + 



M l - i M 2 



U+ + U 



-Mi - i M 2 



T = UTk [-N 2 - i iVi] U + + UTk [-N 2 + i N x ] U. 



+ U 

We notice the properties 



T 



-M 2 + i Mi 



-M 2 + i Mi 



(3.49) 
(3.50) 



S t = E , T T = T . 
The matrix entering in eq. ( p.20| ) can be diagonalized through a unitary matrix C 

CH(k,r])C j = H d (k,ri) diagonal, (3.51) 

such that the energy density is 



H = (a + ,b) H 



a + ,b)H, 



(3.52) 



A first step in this diagonalization can be made by noticing that the matrix H 
can be rewritten as 

H = U ] H U, (3.53) 



with 



hermitean and 



-Mi + iM 2 kN x + ikN 2 
kNt-ikN 2 Mt + iM 2 



U 



( 



-iU* 
-iUl 



(3.54) 



(3.55) 



V2\-U- 

unitary, as it follows from eqs. ( |3.48|) . 

We are not able to provide general analytical formulae for the diagonalization of 
the remaining matrix H . This diagonalization can however be performed numer- 
ically. In addition, some important conclusions can be drawn from the properties 
of the matrix Ti. Due to the relations ( |3.50| ), one can show (i.e. by counting the 
number of independent equations that must be satisfied) that the matrix C entering 
in eq. ( |3.51| ) can be of the form 

I J \ , 

(3.56) 



(3.57) 



( ' iJ* -iP 
(where / and J are N x N matrices). Unitarity of C requires 

J/t + jjt = 1; Pj=J T I*. 
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By explicitly performing the product (|3.51|) one realizes that the eigenvalues of 
the hamiltonian occur in pairs, and that Hd is of the form 

diag(c<Ji, u>2, ■ ■ ■ , ^n, —u)2, • • • , —ojn)- The eigenstates corresponding to each 
couple icUj are to be interpreted as particle and antiparticle states with the same 
energy. Finally, it is possible to show that particles and antiparticles are produced 
in the same amount. Defining the vacuum state to be annihilated by the initial 
annihilation operators <2j and 6j, we have indeed (we remind that in this expression 
we do not sum over i) 

N l ( V ) = (a\a t ) = (b\k) = {JJ ] ) ll . (3.58) 

We assume that no fermionic particles are present at initial time rj. In our formalism, 
this is equivalent to requiring J(fj) = 0. Notice also that the unitarity condition ( [3. 57j ) 
ensures that the Pauli principle is always fulfilled. 



4. One application: non-thermal gravitino production 

In this section we discuss one application of the above formalism, i.e. non-thermal 
gravitino production in a system with two chiral superfields. This model consists of 
two chiral multiplets coupled only gravitationally. The scalar of the first multiplet is 
responsible for driving inflation, while the one of the second breaks supersymmetry in 
the vacuum. The motivations for this analysis, as well for the specific model chosen, 
are explained in the introduction (see also |23| ) . 

This section is divided in five subsections. In the first one we introduce all the 
quantities relevant for this calculation. Although we have exactly followed the con- 



ventions of ref. |19|], the aim of section |4.1| is to provide a practical self contained 



presentation. In section 4.2 we describe the model that we are considering. We also 



discuss there the evolution of the scalar fields, which constitute the external back- 



ground for the fermionic fields. In section |4.3| we show how to apply the formalism 



of the previous section to the calculation of the abundances of the fermions of the 
theory. Our results are presented in the two remaining subsections. In section |4.4j we 



present analytical results in the case in which supersymmetry is actually unbroken in 
the vacuum of the theory. We show that in this case, gravitinos are only gravitation- 
ally (hence negligibly) produced. Already this consideration suggests that in the class 
of models we are considering (i.e. with the two sectors coupled only gravitationally) 
non-thermal gravitino production may be very inefficient in the realistic situation in 
which the observable supersymmetry breaking (TeV scale) is much smaller than the 
scale of inflation (10 13 GeV). This is confirmed by the numerical results presented 



in section \j.5[ which show that gravitino production indeed decreases as the size of 



supersymmetry breakdown becomes smaller. 
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4.1 Definitions 



We write here the relevant equations of motion for the gravitino field and the 



fermionic particles to which it is coupled. We follow the conventions of ref. [[IS 
The starting action is the one of D = 4, N = 1 supergravity, with four fermion in- 
teractions omitted. For simplicity, we do not consider any gauge multiplet, but for 
the moment we allow chiral superfields to be complex. The lagrangian reads 

er l C = ~M*R- gj (d, 0*) (<9%) — V- 

- X - BP + X - m ipfiR ^ ipvR + \ m* ^ L -f tfj uL - 

- 9i j [Xj $ X* + tiP Xj] - m ij Xi Xj ~ mj t X j + 

+ (2 9j % R 7"" X ] + $ R ■ 1VL + h.c.) . (4.1) 



The first line of eq. (|4.1| ) concerns the scalar fields. The first term is the standard 
one of Einstein gravity, with Mp denoting the reduced Planck mass (Mp ~ 2.4 • 
10 18 GeV) and R the Ricci scalar. Conformal time 77 is used and the Minkowski 
metric is taken with signature — h ++. More explicitly, the metric and the vierbein 
are given by = a 2 ^)?^, = a(r})5^, where a is the scale factor of the Universe. 
We then have some chiral complex multiplets formed by (0j, Xi) an d their conjugate 
((f) 1 , x 1 )- K is worth emphasizing that Xi is a l e ft handed field, while x % a right 
handed one. The left and right projections are Pl = (1 + 75) /2, Pr = (1 — 75) /2. 
The gamma matrices in curved space 7 are related to the ones in flat space 7 by the 
relation 7^ = a -1 7 M , and the realization of the latter that we are using is given in 



eq. (|3.29| ). The Kahler metric is the second derivative of the Kahler potential 



while the scalar potential V is defined below. 

In the second line of eq. ( |4.1| ) we have the kinetic and the mass term for the 
gravitino field. The first one is defined to be 

BP = e- 1 e!»""m v D p $ a , (4.3) 

where the covariant derivative 



1 



D^ v = I K + -0^7™. J 5" u - rj v j ip\ , (4.4) 

contains the spin connection u™ n and the connection T^ u (■y mn = [7 m ,7 n ]/2). The 
mass parameter m is instead given by 

K 



m = e 2M pW } (4.5) 
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and it is related to the gravitino mass by 



m & = \m\ Mp . (4.6) 

We then find the kinetic and mass term for the chiral fermions. The first is 
given by 

D^ X i= ^ + ^r7«*)x« + ^[^^-^^i]x* + l1*^^*- ( 4 - 7 ) 

where T{ k = g~ 1 i 1 g k is the Kahler connection. For what concerns instead the 
fermion masses, we have 

m l = D l m = d l m + "j" 2 d l Km , 

Zi IVl p 

m ij _ ^ D J m = [ 9* + -4to ^'^ J m J - rl J m k . (4.8) 
\ 2M| ) 

We can now write, in compact notation, the scalar potential 

V = -3M P 2 \m\ 2 + m i g~ 1 j i m j . (4.9) 



The last line of eq. ( f4.1|) describes the interactions of the gravitino with the chiral 
fields (i.e. with matter). The field vl is defined to be 

VL = m i Xi + {^ i ) X ig.\ (4.10) 



As understood in refs. |18|, [B| ; this combination of matter fields is the goldstino 
(actually its left-handed component) in a cosmological context, where supersymmetry 
is broken both by the kinetic and the potential energies of the scalar fields. We work 
in the unitary gauge, where the goldstino is gauged away to zero. We also Fourier 
transform the fermion fields in the spatial direction, i.e. x{Vi%) = x{v) elx %ki - 

The gravitino field has transversal and longitudinal components. To appreciate 
their different behavior, one can introduce the projectors 



<*-), = 5 (f-ljM*^) 



2 k z 

where k{ are the spatial components of the comoving momentum of the gravitino 
(i.e. d ki = 0) and k 2 = k{ fcj. These projectors are employed in the decomposition 

^ = ipj 1 + (p 7 ) i e + (P k )i h ^ , (4.12) 
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where 9 = 7* ipi . 17 As it is shown in ref . [ 19| , from the lagrangian (|4.1j) one recovers 
four independent equations for the gravitino components. Two of them are algebraic 
constraints which involve ipo, kiipi, and 9. We use them to eliminate the first two 
combinations in favor of the last one. The other two are instead dynamical, and can 
be written in the form 



no . j , a 7 am 



^ T = 



do + B + if hfA )6 rT = 

a a 



(4.14) 
(4.15) 



where 



T = (Xi #o <^ + X J #o <f>i) 
m = Pr m + Pl m* , 



A — - (ax- 7°a 2 ) 



-Ml (3H 2 + 2H 



ni 



B 



3 1 

— aA-\ 7 am7° ( 1 + 3 A 

2 2 M I — 1 1 



«2 = 2 a 1 9q W 



(4.16) 



/ = a _1 9o/, and where H = a/a is the Hubble expansion rate. 

We notice that the transverse component of the gravitino, ipj , is decoupled from 
the longitudinal component and from matter, apart from gravitational effects due to 
the expanding background. In particular, transverse gravitinos are produced only 
gravitationally |14| , 15, [16], and for this reason we will not consider this component 
any longer in the rest of the work. 

We are thus left with the gravitino longitudinal component, rewritten in terms 
of 9, and the matter fields. In case of only one chiral supermultiplet the combination 
T defined above is proportional to the goldstino, and thus vanishes in the unitarity 
gauge. In the more general case of iV chiral superfields, we have (always in the 
unitary gauge) N — 1 non vanishing independent fermionic chiral fields, and one 
should go into a basis orthogonal to the goldstino. The equations of motion for all 
these fields can be of course deduced from the initial lagrangian (|4.1| ). If only two 
superfields are present, one is just left with the matter field T defined above. 



17 Noticc that 



hiP^^O, f(P 7 ) t = l, 
h(P k ) t = l, f(P fc ).=0, 



(4.13) 
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4.2 Description of the model considered and evolution 
of the scalar fields 



The matter content of the model we are considering is of two superfields <3? and S 1 , 
with superpotential 

W = ^$ 2 + yU 2 (/3 + S) (4.17) 
and minimal Kahler potential 

K = &$ + S*S. (4.18) 

The potential for the scalar components and s of the superfields <3> and S can 
be computed using eq. (|4.9| ). We then assume that the scalar fields are real, that is 
(after V is computed) we perform the substitutions 

* = *'^T2' S = S '^7T (4 ' 19) 

In this way the real scalar fields have canonical kinetic terms. 

During inflation, the field acts as the inflaton, while the v.e.v. of s is quickly 
driven to (s) ~ 0. The potential is then practically the one of chaotic inflation, 
and ~ 10 13 GeV must be posed to match the COBE results for the size of the 
CMBR fluctuations. 18 

At the end of inflation, the field oscillates about the minimum = 0. The 
amplitude of these oscillations is dumped by the expansion of the Universe (and, later 
on, also by the decay of the inflaton that every realistic model must include) . If only 
was present, we eventually would have unbroken supersymmetry in the minimum 
= 0. The role of the s field is to provide the supersymmetry breaking in the 
minimum. The second term in eq. ( 4.17Q is known as the Polonyi superpotential [p7 



and it provides a simple example on how supersymmetry can be broken in a hidden 
sector and transmitted to the observable one by gravity. What is remarkable with 
this potential is that, for particular values of the parameter (3, supersymmetry is 
broken with a vanishing value for the cosmological constant. If indeed we take 
(3 = (2 - V3)M P (for a more detailed discussion, see for example ||28||), the potential 
V(4> = 0, s) vanishes in its minimum at 

s = V2(V3 - 1)M P . (4.20) 



18 As it is known, the contributions from the Kahler potential to the scalar potential are very 
relevant for (f> ~ Mp. This is a common problem for supersymmetric theories of inflation, where 
the F terms generically spoil the flatness of the potential during the inflationary regime (for a 
review, see |p5[ ; sec also |^6| for a recent discussion). As a consequence, the theory that we are 
here describing should be modified during inflation; however, we will not consider this issue here 
and we will still assume that the value for is not too different from the one imposed in "usual" 
chaotic inflation. 
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Figure 1: Evolution of the two scalar fields cfr and s for fi 2 = 10 
For this value, the gravitino mass m = e K / 2M ? w evaluates to 



,2 ,,2 



rriQ = e 



^TT-^TT* (4.21) 



M P M P 



which is a typical result for this breaking of supersymmetry. We see that the "inter- 
mediate" scale \x must be taken of order 10 10 GeV to reproduce the expected gravitino 
mass ttLq ~ 100 GeV. 

In the following, we discuss in more details the evolution of the two scalar fields. 
To do this, we use physical time t and work with the adimensional quantities 

i = ±_ % -±_ Q = J_ 

° .\/p ' S Mp ' M v ' 

fstm " a ^a^- ^ s ^- ^ s sfe- (4 ' 22) 

where we remind that H and V are, respectively, the Hubble constant and the scalar 
potential. In terms of these redefined quantities, the equations of motion for the two 
scalars read 

d 2 4>j - ddj dV - . 

+ 3 H — ? + — - = , fa = (f>, s. (4.23) 



dt 2 dt 

We start our numerical calculations at ~ 1.4, short after inflation, and with 
the scale factor a normalized to one. 
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We show in figure [I] the evolution for the two scalar fields after inflation, in the 
case fi 2 = 1CT 2 . As we said, initially the model reproduces the scalar potential of 
chaotic inflation, and thus we have 

S ~ . (4.24) 

The initial dynamics of the Polonyi field s is determined by <fi. More precisely, we 
can write an effective potential V(s) for it by substituting eq. ( |4.24| ) into V(0, s) and 
then averaging over the inflaton oscillations. Expanding V for both and s smaller 
than one, we find that the potential is minimized by 

g ^A 2 [16/^ 2 -(0 4 )] ^ 3^2/3 (ft 2 *) 2 
4(0 2 ) - 16/5 2 A 4 1-Sp 2 (fi 2 t)' 

To be precise, the Polonyi field is always smaller than s mm , due to the fact that the 
expansion of the Universe slows its motion towards the minimum of V(s). However, 
eq. ( 4.25| ) gives a good estimate for the order of magnitude of s in this initial stage. 



What is most important to emphasize, is that eqs. ( 4.24 ) and ( 4.25|) explicitly 



show the presence of two very different (physical) time-scales in the model we are 
considering. The first of them is set by the inverse inflaton mass m^ 1 which is the 
time-scale of the oscillations of the inflaton field. The second one is given by ft 2 m^ 1 . 
Equation ( 4.25|) shows that this is the relevant time-scale for the Polonyi field in the 



initial stage. However, this is true also for the complete evolution of s. To see this, 
let us consider the latest times shown in figure [1|. In this stage the amplitude of the 
oscillations of (f> are negligible. The evolution of s is not any longer influenced by the 
inflaton field, but it starts oscillating about the minimum of its own potential given 
in eq. ( pl20|) . 19 The amplitude of these oscillations is also dumped by the expansion 



of the Universe, while their period is related to the inverse Polonyi mass, which is 
now (i.e. at = 0) given by pF 



m K 



2Vlm Ci ~2Afi 2 m 4> . (4.26) 



The quantity fi 2 defines the ratio between the two scales. In figure [1| we have cho- 
sen, for illustrative purposes, fi 2 = 1CT 2 . However, this value is unphysical, since it 
would correspond to a too high supersymmetry breaking scale. Indeed, as eq. ( |4.21| ) 
shows, we must require fi 2 ~ 10" 11 , if supersymmetry is supposed to solve the 
hierarchy problem. 

While the size of fi 2 controls the supersymmetry breaking in the vacuum of the 
theory, both scalar fields contribute to break supersymmetry during their evolution. 
In particular, both their kinetic and potential energies contribute to the breaking, 

19 There is of course a possible moduli problem associated with these oscillations. However, we 
do not consider this issue here. 



25 




500 



Figure 2: Relative contribution of the two scalar fields eft and s to the supersymmetry 
breaking during their evolution. As in figure [l], fi 2 = 10 -2 . 



as emphasized in ref. |[L8| . This can be seen by looking to the transformation law 
of the chiral fermions Xi under an infinitesimal supersymmetry transformation with 
parameter e. In our case they read 

I 



rrii 



_o "Vi 

71 7 ~dt 



(4.27) 



where fa = <f>, fa = s. 
Following refs. 



19fl ) we define the quantities 

WHS)'- (4 - 28) 

which give a "measure" of the size of the supersymmetry breaking provided by the 
F term associated with the z-th scalar field. More precisely, we will be interested in 
the normalized quantities 



n 



r s = 



(4.29) 



* ft + fi ft + fi 

which indicate the relative contribution of the two scalar fields and s. 

In figure |2| we have shown the evolution of and r s for the specific case 
fi 2 = I0 -2 . As expected, in the initial stages only the inflaton contributes to the 
supersymmetry breaking, while only the Polonyi contributes at later times. The 
regime of equal contribution is around t = fi~ 2 , when <\> and s are of the same size 
(cf. figure |l|). As it should be clear from the above discussion, and r s share the 
identical behavior for all the choices of ft 2 , once t is given in units of pr 2 . 
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4.3 Effective fermionic lagrangian and hamiltonian in the case 
of two chiral supermultiplets 

The fermionic content of the model we are considering is of the gravitino ip^ and 
the two chiral fermions and s. In the unitary gauge, one combination of <ft and 
s, the goldstino v, is set to zero, while the transverse component of the gravitino, 
ipf , is only gravitationally coupled to the other fields. The other two fermions 9 
(the longitudinal gravitino component) and T (the combination of chiral fermions 
orthogonal to v) are coupled together, as we described in section fO|. 

With some algebra, we can rewrite the initial lagrangian ( |4.1| ) in three terms 



C 



^background ^-"ifiT &6 T 



(4.30) 



The first term governs the dynamics of the scalar fields and of the scale factor of 
the Universe. The second describes the (decoupled) transverse gravitino component, 
while the third one reads 



T 



a 



4 k 2 



a 3 9 



fd 9 + if kiA9 + 



- -a7 u + 7 

1 2 ' 2Ml 



am) A6 



am 4k _ 
9 



2 Ml 



a a 



7 U T 



4a 
o^A2 



T 



+ 



fd Q T-if kiAT- -d7°iT + ^AamT + 

10 1 2 ' 2 Ml 



(4.31) 



We have expanded the fermions into plane waves Xi(m k) = e ' iX Xi{r]) } where fcj is 
the comoving momentum (i.e. do ki = 0), and we have introduced 



A = 



a 
2 

a 



<pi fa m k mi (g 1 kl g {j - 5 ik %) 



1/2 



m X (p 2 - m 2 <pi , 



(4.32) 



where the second equality holds in the case of a minimal Kahler potential, g\ = 5j. 
The quantity A has no counterpart in the one chiral superfield case, and indeed it 
is negligible unless both the scalar fields give a sizeable contribution to the breaking 
of supersymmetry. 

One can explicitly verify that the lagrangian ( |4.31]) reproduces the equation of 
motion ( 4.15|) for the longitudinal gravitino component, as well as the one for T 



that one obtains from the initial lagrangian (|4.1|) . However, we notice that the 
two fields 9 and T are not canonically normalized. Canonical normalization has to 
be imposed, if we want our fields to give invariant quantities (as for example the 
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occupation number) in comoving units in the adiabatic regime. Among the possible 
redefinitions, we choose 



(acfif 2 

T = ^(-) T, (4.33) 



2 \aJ 

since the equations of motion look quite symmetric in terms of the new fields. In 
matrix form, they are exactly of the form (|3.33|) , i.e. 

(7 d + i f kiN + M) X = , (4.34) 

where X is the vector (9, Y) T . In our specific case, the "mass" matrix M is given by 

f ma 3 /mo . 

M = diag fe + 2fe ai + aa2 

ma 3 / ma _ . . , , _ , 

+ o T72 «i + « «2 + a (mu + m 22 ) , (4.35) 



2MI 2 \A/| 



and the matrix by 



In the above equations, we have defined = oti/a. The relation a\ + a 2 = 1 
which holds in the one chiral field case [TJ, I3| is now replaced by 20 

a\ + a 2 2 + A 2 = 1 . (4.37) 

We thus see that the matrices N\ and N 2 satisfy both conditions (|3.35|) . 

The equations of motion (|4.34|) have a clear behavior in the low energy limit, 
when the two scalars of the theory settle to their minima. In this final stage one has 
a\ = —1,52 = A = 0, as it can be easily checked from the definitions listed above. 
As a consequence, eqs. ( |4.34j ) decouple, and each of them acquires the standard form 
for spin 1/2 fermions 

(j°d + ifki + am^e = 0, 
(fdo + if h + am f ) f = 0, (4.38) 

where the two masses are constant. In particular, notice that m^ = m/Mp, which is 
exactly the expression that one encounters in supergravity for the gravitino mass. 



20 When only one scalar field gives a substantial contribution to supersymmetry breaking, the 
quantity A almost vanishes, and the relation cti+a^ — l holds approximatively. 
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We thus see that the system has all the properties assumed in section ( |3.2| ), 
so that we can apply the procedure derived there to quantize it and to define the 
occupation numbers of the fermionic eigenstates. Among the possible choices for the 
transformation matrix $ which enters into eq. ( |3.36| ), we take 

$ = -(arccos a x )(^ (4.39) 

with 



u = Jl -a\ = yja% + A 2 . (4.40) 
Following eq. (|3.44 ), the hamiltonian of the system is instead given by 21 

H = x\if kiN + Mi + 7° M 2 ] X , (4.42) 



with 



2 f) ~ 

M x = M+<^-(QMQ-M) + ^Q-^& 1 d Q, 
Z Z UJ z 

M,_ = i^i [M, Q] + doQ] . (4.43) 



A Jl..}- (4-44) 



We have denoted 

A/uj —a 2 /uJ 

As we have already remarked, at late times &x = —1, while a 2 = A = 0. In 
this regime the above hamiltonian becomes the standard one of two decoupled spin 
1/2 fermions 

H = X[ifki + M]X, (4.45) 



with the standard gravitino mass for the field 6 (cf. eq. ([O 

We conclude this subsection discussing the explicit diagonalization of the hamil- 
tonian, i.e. of the matrices H and H entering in eqs. ( |3.52|) and (|3.53|) . One can now 



explicitly verify that the eigenvalues of the Ho matrix occur in pairs, that is they are 
of the form ±u>i, ±u 2 . One can also verify that if (vi, v 2 , v%, is an eigenvector of 
H belonging to the eigenvalue u, then (—v^, —v* A) v *, u|) is also an eigenvector of 



21 Notice that the matrix that appears in the equations of motion (3.38) reads 



M 2 = [Q,M] | + QqL_£1. (4.41) 

At late times M2 ~ Q Q does not vanish. Indeed &2 and A decrease at late times in a way such 
that the elements of Q keep on oscillating with amplitude equal to unity. Therefore, as discussed 



in the footnote before with eq. (3.41), the fields X are not a suitable basis for the definition of 
the hamiltonian. 
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H belonging to the eigenvalue —u. We then find 

rfH R = H d , R=(* ~*}), (4-46) 
diag(c<Ji, UJ2, —0J\, — 0J2) is the matrix that we formally introduced in 



where Hd 
eq. fH). 

The 2x2 matrices defined by eq. ( 3.56|) are thus given by 



V2 



R\ U+ + R\ U. 
-R\ U*_ + R T 2 U 



(4.47) 



We remind that the matrices R\ and R2 are obtained through the diagonalization 
of H see eq. ( |4.46|) . The matrices U + and £/_ are instead determined by their 
evolution equation ( |3.47| ). The only point left is to give more explicitly their values 
at the initial time rj. This can be done by setting J = in eq. ( }4.47| ), which, as we 
remarked, corresponds to requiring no fermions in the initial state (see eq. ( |3.58|) ). 
Moreover, conditions (|3.48|) have to be imposed. From these requirements, we see 



(4.48) 



that U + {rf) has to fulfill 

Ul (fj) \t + R 2 (rj) Rt l (rj) Rj' 1 (rj) R ] 2 (rj) U + (fj) = 2 1 . 



In this last expression, the matrix in square brackets is hermitean and can be diag- 
onalized with a unitary transformation. More precisely, we can set it to be equal to 
WAV with A diagonal and real, and V unitary. The initial condition for U+ can 
thus be written 

U+(rj) = V ] V2A- 1 . (4.49) 
Finally, U_{rf) is obtained by setting J(fj) = in eq. (|4.47| ). 



4.4 Analytical results with unbroken supersymmetry in the vacuum 

The case fi 2 = is particularly interesting since some results can be worked out 
analytically, and since it provides some hints between the final gravitino abundance 
and the size of supersymmetry breaking. For fi 2 = supersymmetry is unbroken in 
the minimum of the theory, at = s = 0. 22 Because of this, in the vacuum of the 
theory the gravitino has only the transverse component. 

The computation of the formulae of section |4.3| is in this case particularly simpli- 
fied. The quantity A vanishes identically, so that the two fields T (which is always 
the Polonyi fermion) and 9 (which is always the inflatino) are decoupled. Going 

22 Indeed for fi 2 strictly zero the potential for the Polonyi held becomes hat for <fi = 0, and s = 
for the whole evolution. 
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back to the formalism of section (|3.2|) , we find that we have to perform only the first 
redefinition of the fermions, X = exp(— 7 $) X, where now $ = diag(<y2, —<f), with 

cos (2 <p) = aci , sin (2 ip) = a 2 ■ (4.50) 

The two redefined fields have the "standard" equations of motion and hamiltonian 

(fd +if ki + m § )§ = 0, 
(fdo+ifh + mf) T = 0, 

H = J d 3 k 9 (i f ki + m § )9 + f (i f k { + m x ) f , 

m § = m e + d (p , rrif = m T — d <p , (4-51) 

with mg and my given in eq. (|4.35|) . 

In practice, "removing" the time dependent matrix which multiplies the momen- 
tum in the original equations for 9 and T gives an additional contribution to the mass 
of the fields 9 and Y, as noticed in |L4L [15]|. When A = 0, an explicit computation 



of <9q0 with the present formalism is provided in |19] 



d (p = -a = a (mn +m)+3a \ H — - m 1 m\ J" 1 — - , (4.52) 
2ai \ V2 J ^(j) 2 + m{ 



where the various quantities have been introduced in section [11] . Using the equation 
of motion for the inflaton field 0, one can show that it is precisely d^ip = mx, so that 
the field T is effectively massless. The mass for 9 is instead of the order the inflaton 
mass. More precisely, it has a variation of the order within the first oscillation 
of the inflaton (i.e. in the time m^ 1 ) and then it stabilizes at fT4] , [T5| . Since 
the fields are decoupled, the formulae for the occupation numbers ( |3.17| ) are quite 
simple. They show that the Polonyi fermion is not produced, while the production 
of the inflatino field has a cut-off at k ~ m^, and decreases as k~ 4 at big momenta. 23 
The main point of this subsection is that the Polonyi fermion is not produced 
at preheating for ft 2 strictly zero. When fi 2 7^ the Polonyi fermion provides the 
longitudinal component for the gravitino, so its abundance turns out crucial to un- 
derstand whether gravitinos are or are not overproduced. If one believes that the 
limit /x 2 — ^ is continuous, the present analysis suggests indeed that the production 
of gravitinos should become smaller as fi 2 decreases. Although we do not have a 
rigorous proof of this continuous behavior, 24 the numerical results that we show in 
the next subsection strongly support this assumption. 



This can be explicitly seen by integrating eq. ( |3.17 ) for B in the limit of large k and with A ~ 1 



The problem is that the dynamics of the Polonyi field is governed by the timescale fi 
which becomes infinite in the limit /t 2 — * [0| . 
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4.5 Numerical results with broken supersymmetry in the vacuum 

We now analyze the situation fi 2 ^ 0. As we have said, in this case the quantity A 
is generally also non vanishing in the most interesting part of the evolution. As a 
consequence, the dynamics of the fermionic fields 9 and T is coupled, i.e. we have 
mixed terms in their equations of motion ( f4.34|) and in their hamiltonian ( }4.42|) . For 
the following discussion it is useful to explicitly write T in terms of the chiral fields 
Xi and X2- Combining the definitions ( |4.16| ), ( }4.32|) , and ( }4.39|) we have, for minimal 
Kahler potential and real scalar fields, 

T = T 9 To 0iXi + 02X2 , (4.53) 

mi 02 — m 2 01 v J 

where we remind that the two scalars 0i and 02 are the inflaton and the Polonyi 
field, while Xi an d X2 the corresponding fermions. Moreover, the definition ([4.1 Of) of 
the goldstino now reads 

v L = - 7 a (p^j Xi ■ (4.54) 

Let us first consider the initial and final stages of the evolution, where only one of 
the two scalar fields significantly contribute to the supersymmetry breaking and the 
quantity A is practically vanishing. During inflation, one has m 2 , 0i, 02 — 0, and 
the supersymmetry breaking is provided almost completely by nil. Moreover, the 
goldstino is practically the field xi- We remind that we are working in the unitary 
gauge, so that Xi — v — 0- Equation ( |4.53| ) thus rewrites 



T- a3/2| 7 ll Xi02 = « 3/2 Xi, £«/T 2 . (4.55) 
mi 02 

Notice the factor a 3//2 appearing in the last expression, which is a consequence of the 
fact that the field T is canonically normalized in comoving units (cf. the discussion 
after the lagrangian ( f4.31|) ). 25 In the late stages of the evolution, supersymmetry 
is instead broken by the Polonyi field, and the only non-vanishing contribution is 
provided by m 2 . With the same arguments used to get eq. (f4.55|) , one can show that 26 



V (XX2 



0, f = -a 3/2 X i, t>/T 2 . (4.56) 



In order to compute the evolution of the occupation number of the fermions we 
adopt the procedure described in section FO: during the evolution of the system, the 



states are mixed in such a way that the hamiltonian is kept in a diagonal form. The 



25 0ne may also worry about dividing by ct>2 — in eq. ( 4.55 ). However, this is due to the fact 
that the quantity T is ill — defined in the static <fii — > limit, while T is not. 
26 0ne can also show that, in the late stages of the evolution, \<j>i\ ~ j^l- 



32 




0.8 f 

0.6 - 

0.4 - 

0.2 - 

^ 



100 200 300 400 500 600 700 800 900 1000 

t 

Figure 3: Evolution of the masses of the two fermionic eigenstates. As in figure [I], 
fi? = 1CP 2 . Notice the different normalizations for the two masses. 



two eigenstates obtained through this diagonalization coincide with the fields 9 and 
T only for A = 0. In particular, this is true at very early and late times. The safest 
way to make the proper identifications in these regimes is to consider the evolution 
of the mass eigenvalues, which always behave like in the example shown in figure ||[ 
The two masses present (for fi~ 2 <C 1) a strong hierarchy. We denote with ip\ the 
eigenstate with bigger mass, and with ip 2 the other one. The mass of ipi converges 
to the inflatino mass (~ 1.31 m^) at late times, and it is always of the order m^. 
On the contrary, the mass of ip2 converges to the gravitino mass (~ 1.31 fi 2 m^) in 
the vacuum. As it will be clear below, we can "qualitatively" identify ip\ with the 
inflatino and ip 2 with the Polonyi fermion for the whole evolution. Although rigorous 
only at late times, this identification can be useful for a qualitative understanding of 
the system. 

What is most important to us is the relation between the eigenstates (ipi, ipz) and 
the gravitino 9 and the matter field T. As we have said, the last fields coincide with 
the physical eigenstates only at the very beginning and at the end of the evolution. 
More precisely, from the behavior of the two masses we have 

B = i) 2 and T = -0i (4.57) 

at late times. On the contrary, it must be 

9 = tpx and T = tp 2 (4.58) 

at early times, since the longitudinal gravitino component 9 is provided by the gold- 
stino and supersymmetry is initially broken only by the inflaton field. 
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Figure 4: Evolution of Nq and Ny for /t 2 = 10 2 and A; = m,A. See the text for details. 

At intermediate times, the hamiltonian cannot be diagonalized with a simple 
rotation in "flavor" space, and 9 cannot be just a simple (i.e. with only numbers 
as coefficients) linear combination of ip\ and ip2- However, we can gain an intuitive 
description of the system through the identifications 



(4.59) 



The coefficients and r s give a "measure" of the relative contribution to supersym- 



metry breaking provided by the two scalar fields (see eq. (|4.29 )). These relations 
can thus be justified as a "generalization" of the equivalence theorem, in a way also 
suggested in [|Tj], [l!J . We remark that they are rigorous at early and late times (when 
they coincide with the identifications ( f4.57|) and ( |4.58 )). At intermediate times they 
interpolate between these two regimes and can be thus used as a qualitative descrip- 
tion of the system. 

From eqs. (f4.59p we deduce the following estimates for the occupation numbers 



N e = r 4> N l + r s N 2 , 
N r = r s N 1 + r N 2 . 



(4.60) 



The evolution of these quantities is shown in figure |] for modes of comoving momen- 
tum k = and for fi 2 = 10~ 2 . Notice that (by construction) Nq = Ni at early 
times, while Ne = A^2 at late ones. In these regimes these identifications are rigorous. 
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Figure 5: Spectrum of the state ipi at different times for p? = 10~ 6 . 
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Figure 6: Spectrum of the state ip2 at different times for fi 2 = 1CP 6 . 

In figures | and || we plot instead the spectra of the states ipi an d ^2 in the case 
fi 2 = and at the times t = 10 m^ 1 (that is, after a couple of oscillations of the 
inflaton), t = fi~ 2 , and t = 10 /i~ 2 m.7 . 

It is apparent that most quanta of the state ^1 are produced at the very first 
oscillations of the inflaton field, while quanta of ?/>2 are mainly produced at the times 
when the Polonyi scalar starts oscillating. This supports the qualitative identification 
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Figure 7: Spectrum of inflatinos at late times. 

of "01 with the inflatino and of ip2 with the Polonyi fermion. It is worth noticing 
that, for comoving momenta smaller than m^, the increase of N 2 is not related to a 
"conversion" of quanta of ip\ to ip2- Indeed, the increase in N 2 (k) is not accompanied 
by a decrease in Ni(k) for k < m^. 

We can now show our most important result, that is the spectra of T and 6 
at the end of the process. We present them in figures [7| and respectively. They 
are computed at the time 27 t = lOjl^m^ 1 . The time required for the numeri- 
cal computation increases linearly with fi 2 , and the realistic case fi 2 = 10 -11 is far 
from our available resources. We thus kept (i 2 as a free parameter and we per- 
formed the explicit numerical computation only up to fx 2 = 10~ 6 . In particular, in 
figures ^ and |8| the spectra of the fermions produced at preheating are shown for 
jj 2 = lO" 2 , 1(T 3 , 10 -4 , 1(T 5 , 10 -6 . The case fi 2 = l(T n can be clearly extrapolated 
from the ones shown in these figures. 

In figure |7], the spectra for the state ip\ are shown. This state corresponds to 
the matter fermion T in the true vacuum. It is apparent that the main features 
of the spectrum are independent of the value of fi 2 . The reason for this is that ipi 
is associated to the inflatino, that is produced by the coherent oscillations of the 
inflaton. As we discussed above, the dynamics responsible for the production of this 



27 In the cases fi = 10~ 2 — 10~ 4 we have continued the evolution further, until the spectra stop 
evolving. We have found that the spectra shown in figure coincide with the final ones, while Ng 
very slightly decreases for i > 10 /t — 2 ml 1 . Thus, we believe the results shown in figure ^ to provide 
an accurate upper bound on the final gravitino abundance. 



o 
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O 
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Figure 8: Spectrum of gravitinos at late times. 

state is independent on the value of /t 2 . Therefore the only relevant scale for N\ 
is m^, and indeed the spectrum of ip\ exhibits a cut-off at comoving momentum 
k ~ m^. 

The spectra shown in figure || are related to the abundance of gravitinos after 
the fields have stabilized in their minima. We see that in this case the occupation 
number decreases as fx 2 becomes smaller. Indeed, the occupation number N 2 {k) is 
of order unity for comoving momenta k smaller than some cut-off k*. From figure ||] 
we can deduce the dependence k* oc (/i 2 ) 1//3 of this cut-off on the parameter jl 2 . This 
behavior suggests that gravitinos are not produced in the limit (x 2 — > 0, confirming 
what we have argued in the section P . 

We can conclude that in the model we are considering both inflatinos and grav- 
itinos are produced nonthermally. However, the mechanism responsible for the pro- 
duction of gravitinos is much less efficient that the one acting on inflatinos. Indeed, 
the latter is related to the dynamics of the inflaton, while the former is related to 
the dynamics of the Polonyi field. As a consequence, and as it is clearly confirmed 
by the comparison of figures [7| and |8], the number of non-thermal gravitinos is much 
smaller than the number of inflatinos. 



o 

O 
O 



o 
o 



5. Conclusions 

Particles with gravitational decay constitute one of the most serious danger for cos- 
mology, since they can spoil the successful predictions of primordial nucleosynthesis. 
This is particularly true for gravitinos in models where supersymmetry is gravita- 
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tionally broken. To avoid this danger, strong upper limits must be imposed on their 
number. While gravitino thermal production has been very well studied and under- 
stood over the last twenty years, only very recently non-perturbative creation after 
inflation has been considered. Non-perturbative production appears more involved 
than the thermal one, and indeed explicit calculations were so far performed only 
in the one superfield case with supersymmetry unbroken in the vacuum. The re- 
sults of these calculations showed that the coherent oscillations of the inflaton field 
easily cause an overproduction of the longitudinal gravitino component. However, 
this component is absent in the vacuum of the theory if supersymmetry is unbroken. 
One is thus led to wonder if these quanta are really dangerous gravitinos or should 
be better understood as harmless inflatinos (in case of only one chiral multiplet, 
the scalar is necessary the inflaton). In order to discriminate between inflatino and 
gravitino production it is necessary to consider more realistic schemes. The simplest 
natural possibility is to consider two separate sectors, one of which drives inflation, 
while the second is responsible for supersymmetry breaking today. In the present 
analysis we have studied the situation in which the two sectors communicate only 
gravitationally, in the hope that suppressing their interaction may lead to a small 
gravitino production. 

The presence of more than one chiral superfield increases the difficulty of the 
problem. One has to deal with a (quite involved) coupled system, constituted by the 
gravitino longitudinal component and some combinations of the matter fields, in the 
external background made by the scalar fields of the theory. The most difficult part 
of this analysis is to provide a formalism in which the coupled system is quantized, 
with a clear definition of the occupation numbers for the physical eigenstates. This 
is a very interesting problem by itself, which can have several other applications 
besides the gravitino production. We faced this problem in the first part of the work. 
We showed that the standard procedure for quantizing and providing the occupation 
number for one field can be generalized to systems of multi coupled fields, both in 
the bosonic and in the fermionic case. Although far from trivial, this generalization 
can be presented in a remarkably simple form. The application of this formalism 
to the gravitino production is performed in the last section of the paper. We have 
shown that, in the specific model considered here, the number of produced gravitinos 
is very sensitive to the size of the final supersymmetry breaking, and that it actually 
vanishes in the limit in which supersymmetry is preserved. Due to the small scale 
of the expected supersymmetry breakdown, we conclude that gravitino non-thermal 
production is very suppressed in this model. 

Note added. Contemporarily to the present manuscript, the work |29| appeared 
on the database. This work does not consider gravitinos, but studies the non-thermal 
production of moduli fields coupled to the inflaton sector. Moduli fields can also be 
very dangerous for cosmology because, analogously to gravitinos, they are expected to 
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decay gravitationally. The work |29| shows that moduli production can be significant 
if — after inflation — they are strongly coupled to the inflaton. These considerations 
seem to agree with our suggestion p3| that strong couplings between potentially 
dangerous relics (in our case gravitinos) and the inflaton sector should be avoided in 
any viable cosmological model. 
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